% In the meantime, here is a file of the CHAMP EEJ currents 
% in which I have added the F10.7A as an additional column. 
% F10.7A is the 81 day average of the F10.7. As Tim Fuller-Rowel 
% explained to me, the best proxy for the EUV radiation is NOT F10.7, 
% but 0.5*(F10.7+F10.7A), also sometimes called EUVAC index because 
% it is used in the EUVAC model. 

%solar flux is the  Absolute daily solar radio flux (W/m^2/Hz) obtained
%from spidr website

load  C:\Manoj\projects\plasma\solar_flux.mat;

f107a = runmean(solar_flux_new,40,[],'edge');
euvac = 0.5*(solar_flux_new+f107a);

save C:\Manoj\projects\ace\solar_flux_new.mat solar_flux fday f107a euvac;


% The foolowing lines padd 8 days to euvac, find common days with Julia
% observations and save
solar_flux(end+1:2990) = 61.7;
euvac(end+1:2990) = 61.8870;
solar_fday = [solar_fday [ 2982:2989]];
[c,ia,ib] = intersect(solar_fday,floor(fday(:,1)));
julia_solar_flux = solar_flux(ia);
julia_euvac = euvac(ia);

